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Active nanoplasmonic metamaterials support bright and dark modes that compete for gain. Using 
a Maxwell-Bloch approach incorporating Langevin noise we study the lasing dynamics in an active 
nano-fishnet structure. We report that lasing of the bright negative-index mode is possible if the 
higher- Q dark mode is discriminated by gain, spatially or spectrally. The nonlinear competition 
during the transient phase is followed by steady-state emission where bright and dark modes can 
coexist. We analyze the influence of pump intensity and polarization and explore methods for mode 
control. 
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Plasmonic nanolasers, 1 4 arguably the smallest attain- 
able laser devices, are precursors towards active inte- 
grated nano-optics. Owing to their ability to bind and co- 
herently amplify light, these subwavelength devices may 
also be regarded as coherent "meta-atomic" sources, el- 
ementary building blocks of light-emitting metamateri- 
als. In isolation, a nanolaser utilizes feedback by stim- 
ulated emission into a strongly localized, weakly ra- 
diative resonance of a metal-dielectric nanoparticle 2 or 
nanostructure. 3,5 For deep-subwavelength particles ra- 
diative loss becomes negligible and the nanolaser res- 
onates with modes that are effectively trapped on the 
metal surface. In this limit, the nanolaser is a spaser, 6 
a dark coherent emitter of plasmons that can be theo- 
retically described in quasistatic approximation. Conse- 
quently, the spasing threshold is reached when gain over- 
comes dissipative loss and only depends on the intrinsic 
dispersion of the materials but not on the geometry. 7 ' 8 

The introduction of gain into optical metamaterials, on 
the other hand, has been primarily motivated by the aim 
to design loss- free metamaterials. 9 13 Beyond the point of 
loss compensation, when gain outweighs dissipative loss, 
the structure becomes amplifying. 14 Higher gain densities 
(and pump intensities) may lead into the lasing regime if 
dark modes are suppressed and dissipative and radiative 
losses are overcome. 15,16 The required feedback emerges 
from the coupling of the bright and dark meta-atomic res- 
onances to Bloch bands cj^(fc), which can differ greatly 
in terms of their damping constants. With increasing 
pump strength, and depending on the band dispersion 
and the gain spectrum g(w), an increasing number of k 
states becomes unstable (freely oscillating), entering the 
nonlinear competition for gain. If dispersion is weak (flat 
bands), the k states will equally benefit from the gain 
and spatial hole burning effects will dominate, causing a 
complex spatio-temporal light-field dynamics known as 
filament at ion. 17 Hence, two aspects require clarification. 
First, more fundamentally, whether coherent light emis- 
sion is possible in nanoplasmonic gain-enhanced metama- 
terials despite the unique presence of dark modes, which 
may clamp the gain before any bright mode can cross 
the threshold. Beyond that a second, technically rele- 




FIG. 1. (color online) Electric field enhancements and charge 
distributions on the metal films inside the double fishnet unit 
cell for the bright negative-index mode (left) and the dark 
quadrupole mode (right). The direction of the electric field 
in the x-z plane is indicated by the white arrows. 



vant, aspect concerns the control of spatial coherence by 
engineering the dispersion of the uj^(k) bands to effec- 
tively suppress certain Bloch states, enabling, e.g. narrow 
angle emission from large-area surfaces. 

In this Rapid Communication, we investigate the fun- 
damental nonlinear dynamic interplay of bright and dark 
lasing states in an active nanoplasmonic metamaterial 
(see Fig. 1). To obtain insight into the nature of this com- 
petition, we reduce the potentially overwhelming com- 
plexity of the dynamic many- mode interplay to fcy = 
(normal to plane) emission, where the dipolar mode is 
brightest and the quadrupolar mode is truly dark (does 
not radiate). A four- level Maxwell-Bloch approach for 
gain-enhanced nanoplasmonic metamaterials, previously 
employed to study loss-compensation in nano-fishnet 
metamaterials in an ultrafast pump-probe setup 10 ' 18 and 
in an amplification setting, 14 is extended by stochastic 
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Langevin forces 19 23 to accommodate for bath-induced 
quantum noise. 24 This microscopic ab initio approach 
takes into account all relevant scattering channels and 
incorporates the two prevalent forms of light and plas- 
mon emission, amplified spontaneous emission (ASE) and 
coherent stimulated emission, into a self-consistent theo- 
retical framework. 

Our investigations focus on a double-fishnet structure 
with unit-cell side length of p = 280 nm and hole side 
lengths of 120 nm and 100 nm in x and y directions, re- 
spectively. The thickness of the two holey silver films 
is 40 nm and the spacer layer between the films is 60 nm 
thick. Both the spacer layer and the holes are filled with a 
dielectric of refractive index n = 1.62. The (cold-cavity) 
fey = modes of the passive structure are a bright mode 
at 713.8 nm with Q factor Q = 52 and a dark mode 
at 732.2 nm with Q = 134. Figure 1 depicts the en- 
hancement and direction (white arrows) of the electric 
field for both modes in the x-z plane. A further analysis 
of the passive modes can be found in the Supplemen- 
tal Material. 20 ' 25 27 Whilst the two modes have different 
shapes, they overlap in the vicinity of the hole and will, in 
active structures, compete for gain. To explore the non- 
linear dark-bright mode dynamics above threshold, gain 
is introduced into the spacer region and the holes of the 
structure. We describe the gain medium by a four-level 
system absorbing at 680 nm with a coupling strength of 
a a — 1.3 C 2 kg _1 and emitting at A e = 718 nm with a cou- 
pling strength of a e = 1.03 C 2 kg -1 . The density of the 
gain molecules is set to N = 3 TO 19 cm -3 and the spectral 
width of the emission line is chosen to be Y e = 1/ (20 fs); 
further gain parameters and a detailed description of the 
four- level model have been presented elsewhere. 18,20 We 
assume that the system is continuously pumped with 
^-polarized (later ^-polarized) light at a wavelength 
of 680 nm and that at t = all gain molecules reside 
in their ground state. In a first numerical experiment, 
we study the transient behavior at a pump intensity of 
10.75 MWcm -2 , well above threshold. Figure 2 shows the 
temporal evolution over the first 90 ps of the total signals 
(black areas), the filtered signals (solid white, dashed 
green and dashed-dotted yellow lines) in both far-field 
[Fig. 2(a)] and near-field [Fig. 2(b)], as well as the aver- 
age inversion [red solid line in Fig. 2(b)]. At early times 
(t < 17 ps), the total far-field intensity, composed only 
of the transmitted pump field, steadily rises due to satu- 
ration of the absorption. During this phase the internal 
inversion also continuously rises. After approximately 
18 ps, a strong burst is observed in the far-field, initi- 
ating relaxation oscillations of the intensity and inver- 
sion on a time-scale of 3.8 ps. Initially, these oscillations 
seem to lead to stable steady-state emission but are then 
disrupted by large intensity dips starting at t = 45 ps. 
Following these dark oscillations, emission only steadies 
again at around t = 80 ps. 

Using a time-dependent frequency filter, 20,28,29 we de- 
compose the far-field signal into its frequency compo- 
nents. Two dominant frequencies are found: The pump 
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FIG. 2. (color online) (a) Far- field intensity in transmission 
and (b) energy and average inversion (red solid line, right axis) 
inside the nano- fishnet metamaterial over time. The signals, 
time-averaged over 0.4 ps (black), are decomposed into the 
pump mode at 680 nm (green dashed line), the bright mode at 
717.25 nm (white solid line), and the dark mode at 731.8 nm 
(yellow dash-dotted line), (c) Inversion profiles in the x-y 
plane centered between the metal films at the times indicated 
in (b). 



at 680 nm (green dashed line) and, subject to frequency 
pulling, 30 the bright mode at 717.25 nm (white full line). 
Figure 2(a) clearly shows that the transmitted pump 
intensity stabilizes as the bright relaxation oscillations 
dampen in the far-field. The origin of the intensity os- 
cillations at 50 ps is not obvious from the far-field mea- 
surements alone, yet we observe that the bright mode 
settles at a lower intensity value after these dips. To un- 
derstand this behavior we resolve the near-field dynamics 
by recording and filtering the temporal evolution of the 
internal energy density [Fig. 2(b)]. This reveals that the 
dips in the far-field are associated with relaxation os- 
cillations of the internal dark mode at 731.8 nm (yellow 
dash-dotted line) with a frequency of ~(15ps)~ . Strong 
mode competition causes the bright mode intensity to 
drop when the dark mode peaks during its relaxation 
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oscillation. As the bright and dark modes undergo relax- 
ation oscillations, the distribution of inversion changes 
due to spatial hole burning effects. In Fig. 2(c) we show 
the profiles in the x-y plane extracted at three different 
points in time: (I) just before the first spike of the bright 
mode, (II) before the first spike of the dark mode, and 
(III) when the bright and dark mode energies inside the 
resonator have settled. Before lasing sets in, inversion 
is deposited in regions where the pump mode has high 
intensity values [Fig. 2(d)] . When the bright mode be- 
comes unstable, inversion is rapidly depleted in regions 
where its intensity is highest, while the pump keeps pro- 
viding inversion to other regions [Fig. 2(c.II)]. Eventu- 
ally, with the inversion still rising, the dark mode re- 
ceives enough gain to become unstable, further reducing 
the inversion around the hole. After this point, the high- 
est inversion is found within a small square area around 
the center, where stimulated emission is weak as both 
modes have a nodal point in the center of the unit cell 
[Fig. 2(c.III)]. 

The fact that the bright mode starts oscillating first 
indicates that this mode has a lower threshold than the 
dark mode. To further understand the threshold behav- 
ior and the relative strengths of the dark and bright 
modes, we numerically extract the steady-state intra- 
cavity energies for varying pump intensity. We first pump 
the system with high intensity, effectively reducing the 
damping time of the relaxation oscillation that is in- 
versely proportional to the pump intensity. By ramping 
down the pump power with increasingly longer steps of 
70 — 210 ps, we ensure that the modal energies can adi- 
abatically follow without pronounced relaxation bursts 
[inset in Fig. 3(a)]. The so extracted input-output curves 
for the two modes, log- log plotted in Fig. 3(a), show an 
asymptotic linear dependence on the pump power above 
and below threshold. Owing to the presence of ampli- 
fied spontaneous emission, a gradual transition to lasing 
occurs for both the bright and the dark mode. Within 
our modelling approach spontaneous emission enters in 
the form of Langevin noise originating from the strong 
system-bath coupling in the gain medium. Indeed, our 
results confirm that the induced dipole fluctuations qual- 
itatively reproduce the threshold behavior predicted by 
(class-B) rate equation models that incorporate sponta- 
neous emission in the form of phenomenological terms, 
which account for the coupling of inversion to the ground- 
state of the quantized photon field. The threshold of 
the bright mode is reached at a pump energy density of 
approximately 0.03 mJcm -3 (corresponding to a pump 
intensity of 0.3 MWcm -2 ). In contrast, the threshold 
of the dark mode is crossed at a pump energy density 
of around O.lmJcm -3 (1.1 MWcm -2 ), more than three 
times higher than that of the bright mode. 

Thus, despite the higher Q factor of the dark mode, 
it is the bright mode that reaches threshold first. This 
implies that, for the chosen set of optogeometric parame- 
ters, the gain provided to the bright mode is significantly 
higher compared to that available to the dark mode. 
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FIG. 3. (color online) (a) Dependence of the average en- 
ergy densities of the bright (red dots) and dark mode (black 
squares) on the pump mode average energy density. The in- 
set illustrates how the energy densities follow the ramp-down 
of the pump starting at the pump intensity marked by the 
arrow, (b) Mode intensity (left), inversion (middle) and local 
gain [right, see Eq. (1) in the x-y plane centered between the 
metal films. 



In general, the local gain is both time- and frequency- 
dependent and can be defined as 



gioc(w,r,t) = 



a e C(uj)\E(u;,r)\ 2 AN(r,t) 
2T e J v u(u,r)d 3 r 



(1) 



where C(uo) is the Lorentzian emission line and u is the 
energy density. Accordingly, the overlap of the intensity 
profile with the inversion strongly affects the gain deliv- 
ered to the modes. In the left column of Fig. 3(b) we de- 
pict the normalized intensity profiles \E(lj, r)| 2 / \E max \ 2 
of the bright (top) and dark (bottom) modes in the 
x-y plane. The plot in the middle shows the steady- 
state distribution of the inversion A7V(r, t —> oo) in the 
same plane. Multiplying the two profiles with the emis- 
sion line results in the local gain profile at steady state 
gioc{u, r, t — » oo) [right column in Fig. 3(b)]. In this way, 
it can be clearly seen that the gain rate, given by the vol- 
ume integral over the local gain, is significantly higher for 
the bright mode, as expected. 




Effective strategies for mode control can be devised by 
closer inspection of Eq. (1), which states that the modal 
gain depends on three factors: the gain spectrum, the 
spatial mode shape, and the dynamic inversion profile. 
By tuning the overlap between the gain spectrum and 
the resonance lines of the bright and dark modes one 
can change their relative intensities. Figure 4(a) depicts 
the gain spectrum (shaded area) together with the reso- 
nances of the bright and dark modes. The initial choice of 
A e = 718 nm for the gain maximum implies a strong over- 
lap of the gain spectrum with the bright mode (713.8 nm) 
and hence a higher gain for this mode. Favored spec- 
trally, the effective gain of the bright mode exceeds that 
of the dark mode, explaining the lower threshold (Fig. 3) 
and why bright emission sets in first (Fig. 2). To exam- 
ine the impact of spectral alignment, we shift the gain 
peak relative to the resonances. Figure 4(b.I) shows the 
variation of the steady-state output of the bright (circles) 
and dark modes (squares) depending on the spectral po- 
sition of the gain peak. We observe a regime of coexis- 
tence of the modes between 713 nm and 723 nm (marked 
by dotted lines). At 723 nm, halfway between the reso- 
nances, the bright mode switches off. Conversely, below 
713 nm, with the gain peak shifted towards the lower- 
wavelength side of the bright resonance, the dark mode 
switches off and only bright emission occurs. Hence a 
strong spectral discrimination is required to completely 
suppress the dark mode owing to its lower losses. These 
results underpin the necessity of spectral engineering for 
mode control, which considers both the structural reso- 
nances and the fluorescence spectrum. In practice, how- 
ever, the gain spectrum is often predetermined by the 
choice of material or application. Instead, one can shift 
the metamaterial resonances by an appropriate choice of 
geometric parameters to achieve spectral mode control. 

Spectral tuning is not the only way of controlling the 
state of emission. The effective gain of a mode is, as 
Eq. (1) indicates, also proportional to the overlap inte- 
gral of the mode intensity profile and the inversion pro- 
file J v d 3 r|£(cj,r)| 2 AN(r,i). While the mode profile 
is, to a good approximation, determined by the passive 
metamaterial structure, the inversion profile additionally 
depends on the pump excitation, in particular on the 
pump's modal field enhancement that strongly imprints 
on the inversion. For example, for the specific structure 
under investigation, the spatial deposition of inversion 
can be controlled by switching from E x - to ^-polarized 
pump light. Figures 4(c.I) and (c.II) depict the small- 
signal inversion profiles in the x-y plane (before lasing 
sets in) generated by E x - and ^-polarized pump light, 
respectively. When the structure is pumped with E x - 
polarized light, the inversion is distributed in and around 
the holes with a profile that closely resembles the butter- 
fly shape of the dipolar bright mode [Fig. 4(d)] . For 
an £^ -polarized pump, on the other hand, the inver- 
sion concentrates mainly in the projection of the holes 
as the pump couples to the extraordinary transmission 
resonance (EOT) at 707.7 nm for this polarization 26 ' 31,32 
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FIG. 4. (color online) Dependence of steady-state emission 
into dark and bright modes on emission wavelength and po- 
larization, (a) The emission line (here displayed for A e = 
718 nm) and dark and bright mode lines are shown schemat- 
ically indicating the maximum emission at A e and how it is 
shifted for control, (b) Average mode energy densities dur- 
ing steady state in dependence of the emission wavelength of 
the dye for E x (I) and E y polarization (II) of the pump field 
with intensity 10.75 MWcm -2 . The bright mode is given by 
the red line and dots, the dark mode by the black line and 
squares, (c) Inversion profiles before the onset of relaxation 
oscillations for E x (I) and E y polarization (II) highlighted by 
the black arrows. 



[Fig. 4(c.II)]. The fact that the bright mode benefits less 
from the inversion deposited by ^-polarized pump-light 
than the dark mode is confirmed by Fig. 4(b.II) where 
we show the steady- state output for £^ -polarized pump- 
ing. A comparison with Fig. 4(b.I) reveals that for E y - 
polarized light the dark mode becomes more intense, the 
bright emission weakens and the window of dark-bright 
coexistence displays a blueshift. 

Finally, we note that the competition between bright 
and dark modes as well as schemes for mode control de- 
pend on the geometry. Using a similar structure with 
square holes, for example, the E x and E y bright modes 
are degenerate and the dark mode is then found on the 
high-frequency side of the bright mode. It is therefore 
possible to completely suppress the excitation of the dark 
mode by overlapping its resonance with the absorption 
line of the gain medium. In this case, dynamic competi- 
tion between the two bright modes of orthogonal polar- 
ization can be observed. 19 

In conclusion, we have analyzed the nonlinear dynamic 
competition of dark and bright lasing states and methods 
of mode control. Using a full-vectorial Maxwell-Bloch 



5 



Langevin approach, strong nonlinear mode competition 
during the transient regime was found, as the modes are 
harnessing gain in disparate, yet overlapping spatial re- 
gions. The steady state input-output curves of the bright 
and dark modes reveal a characteristic threshold behav- 
ior and a regime of stable coexistence. By changing the 
spectral alignment of the emission line or the spatial de- 
position of gain, the relative intensities of the modes can 
be manipulated to the limits where either mode switches 
off, and sole bright or dark emission occurs. Moreover, 
in extending the calculations to multiple unit cells, the 
presented approach will naturally incorporate the com- 
petition between Bloch states, opening the possibility to 



investigate the dynamic build-up of spatial coherence and 
the emerging radiation patterns. Our work may aid the 
design of emitting metamaterials by promoting the the- 
oretical understanding of the dynamic interplay between 
and control of lasing states in active nanoplasmonic struc- 
tures. 
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